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Abstract. We study the binding energy differences of the valence proton and neutron 
of the mirror nuclei, 15 - 15 N, 17 F - 17 0, 39 Ca - 39 K and 41 Sc - 41 Ca, using the quark- 
meson coupling model. The calculation involves nuclear structure and shell effects 
explicitly. It is shown that binding energy differences of a few hundred keV arise from 
the strong interaction, even after subtracting all electromagnetic corrections. The 
origin of these differences may be ascribed to the charge symmetry breaking effects 
set in the strong interaction through the u and d current quark mass difference. In 
this report, we first review the quark-meson coupling model. In particular, we discuss 
about the nucleon mass in nuclear medium. Then, we present details of the charge 
symmetry breaking in finite nuclei, especially the Okamoto-Nolen-Schiffer anomaly. 



INTRODUCTION 

The discrepancy between the calculated binding energy differences of mirror nu- 
clei and those measured is a long-standing problem in nuclear physics. It is known 
as the Okamoto-Nolen-Schiffer (ONS) anomaly [1,2]. Although it was first thought 
that electromagnetic effects could almost account for the observed binding energy 
differences, it is now believed that the ONS anomaly has its origin in charge sym- 
metry breaking (CSB) in the strong interaction [3]. In addition to calculations 
based on charge symmetry violating meson exchange potentials [3-5], a number of 
quark-based calculations have been performed [6,7] in an attempt to resolve this 
anomaly. In such calculations, CSB enters through the up (u) and down (d) current 
quark mass difference in QCD. Despite these efforts, the difficulty of producing a 
realistic description of nuclear structure on the basis of explicit quark degrees of 
freedom has hindered the direct calculation of the binding energy differences. 

In this study we report the results for the binding energy differences of the va- 
lence (excess) proton and neutron of the mirror nuclei, 15 - 15 N, 17 F - 17 0, 39 Ca 
- 39 K and 41 Sc - 41 Ca, calculated using a quark-based model involving explicit nu- 
clear structure and shell effects, namely the quark-meson coupling (QMC) model 
[8] . This model has been successfully applied not only to traditional nuclear prob- 
lems [8] but also to other new areas as well [9]. Although some exploratory QMC 
results on the ONS anomaly have already been reported [7], an early version of the 



model was used there, and it was applied to finite nuclei only through local density 
approximation, rather than a consistent shell model calculation. 

THE QUARK-MESON COUPLING MODEL 

In this section, we introduce the QMC model, and then report the medium 
modification of the nucleon properties in finite nuclei [8]. 

Effect of Nucleon Structure 

Let us suppose that a free nucleon (at the origin) consists of three light (u and 
d) quarks under a (Lorentz scalar) confinement potential, V c . Then, the Dirac 
equation for the quark field ip q is given by 

\i 1 -d-m q -V c (r)]i> q (r) = 0, (1) 

where m q is the bare quark mass. 

Next we consider how Equation (1) is modified when the nucleon is bound in 
static, uniformly distributed (iso-symmetric) nuclear matter. In the QMC model [8] 
it is assumed that each quark feels scalar, V^, and vector, V£, potentials, which are 
generated by the surrounding nucleons, as well as the confinement potential. This 
assumption seems appropriate when the nuclear density p_e is near around normal 
nuclear matter density (p = 0.15 fm" 3 ). If we use the mean-field approximation 
(MFA) for the meson fields, Equation (1) may be rewritten as 

[i 7 • d - (m 9 - V s q ) - V c (r) - = 0. (2) 

The potentials generated by the medium are constants because the matter dis- 
tributes uniformly. As the nucleon is static, the time-derivative operator in the 
Dirac equation can be replaced by the quark energy, — ie q . By analogy with the 
procedure applied to the nucleon in QHD [10], if we introduce the effective quark 
mass by m* — m q — V s g , the Dirac equation (2) can be rewritten in the same form 
as that in free space with the mass m* and the energy e q — instead of m q and e q . 
In other words, the vector interaction has no effect on the nucleon structure except 
for an overall phase in the quark wave function, which gives a shift in the nucleon 
energy. This fact does not depend on how to choose the confinement potential, V c . 
Then, the nucleon energy at rest in the medium is given by En = M N (V£) + 3V£, 
where the effective nucleon mass M N depends on only the scalar potential. 

We can extend this idea to finite nuclei [8]. Let us suppose that the scalar 
and vector potentials in Equation (2) are mediated by the a and uj mesons, and 
introduce their mean-field values, which now depend on position r, by V£(r) = 
g%o~(r) and V£(r) = g^uj(r), respectively, where g q a (g^) is the coupling constant of 
the quark-a (cu) meson. Furthermore, we shall add the isovector vector meson, p, 



and the Coulomb field, A, to describe finite nuclei realistically. Then, the effective 
Lagrangian density for finite nuclei would be given by [8] 

_ T N g 

Cqmc = #7 • d — M% — g^o - 9 P ^h - ^{l + ?f )A 7o ]V> (3) 

- \[(^) 2 + m\A + \[iyuf + ml^\ + \[{vbf + m y] + \iyA)\ 

where ip and b are respectively the nucleon and the p fields. m a , and m p are 
respectively the masses of the a, u and p mesons. g u and g p are respectively the 
u-N and p-N coupling constants, which are given by g w = 3g q and g p = g q p (where 
g q is the quark-p coupling constant). 

If we define the field-dependent cr-N coupling constant, g a (&), by [8] 

M* N (a(f)) = M N - g a (a(r))a(r), (4) 

where is the free nucleon mass, it is easy to compare with QHD [10]. The 
difference between QMC and QHD lies only in the coupling constant g a , which 
depends on the scalar field in QMC while it is constant in QHD. However, this 
difference leads to a lot of favorable results [8]. 

Now we need a model for the structure of the nucleon involved to perform actual 
calculations. We here use the MIT bag model. In the present model, the bag 
constant, B, and the z parameter for the nucleon are fixed to reproduce the free 
nucleon mass (M N = 939 MeV) and the free bag radius R N = 0.8 fm. In the 
following we choose m q = 5 MeV and set m a = 550 MeV, = 783 MeV and m p 
= 770 MeV. (Variations of the quark mass and Rn only lead to numerically small 
changes in the calculated results [8].) We then find that B 1 ^ = 170.0 MeV and z 
= 3.295. 

For infinite nuclear matter, from the Lagrangian density (3), we can easily find 
the total energy per nucleon, E tot /A, and the mean-field values of uo and p (which 
are respectively given by baryon number conservation and the difference in proton 
and neutron densities). The scalar mean-field is given by a self-consistency con- 
dition, dE tot /da = [8]. The coupling constants, g% and are fixed to fit the 
average binding energy (—15.7 MeV) at po for nuclear matter. Furthermore, the 
p-N coupling constant is used to reproduce the bulk symmetry energy, 35 MeV. We 
then find [8]: <? 2 /47r = 5.40, g 2 /47r = 5.31, g\j^K = 6.93, and the nuclear incom- 
pressibility, K ~ 280 MeV. Note that the model gives the variation of the nucleon 
bag radius, 5Rn/Rn = —0.02, the lowest quark eigenvalue, 5x*/x q = —0.16 and the 
root-mean-square radius of the quark wave function, Sr*/r q = +0.02, at saturation 
density. 

Using these parameters, we can solve a finite nuclear system. As an example, 
we show charge density distribution of 40 Ca in Figure 1. The QMC model can 
reproduce the properties of not only nuclear matter but also finite nuclei (for more 
details, see [8]). 



Nucleon Mass in Nuclear Matter 



Here we consider the nucleon mass in matter furthermore. The nucleon mass is 
a function of the scalar field. Because the scalar field is small at low density the 
mass may be expanded in terms of a as 

Since the interaction Hamiltonian between the nucleon and the a field at the quark 
level is given by H int = —3g% J df ip aijj q , the derivative of Af^ with respect to a 
is —3g% J dripgipq = — 3#£5V(cr), where we have defined the quark scalar charge in 
the nucleon, Sn(ct), which is itself a function of a. Because of a negative value of 
the derivative, the nucleon mass decreases in matter at low density. 

Furthermore, we define the scalar-charge ratio, S^ia^j Sn(0), to be Cjv(o") and 
the <7-N coupling constant in free space to be g a (i.e., g a = g a (a = 0) = 3g£Sjv(0)). 
Using these quantities, we find 

M* N = M N - g a a - l - g<7 C' N {0)a 2 + 0(a 3 ). (6) 

In general, CV is a decreasing function because the quark in matter becomes more 
relativistic than in free space. Thus, C' N (0) takes a negative value. If the nucleon 
were structureless would not depend on a. Therefore, only the first two terms in 
the RHS of Equation (6) remain, which is exactly the same as the effective nucleon 
mass in QHD [10]. 
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FIGURE 1. Charge density distribution for 40 Ca [8] compared with the experimental data and 
that of QHD. 



TABLE 1. Inputs, parameters and some of the quantities calculated in the present study. 
The quantities with a star, *, are those quantities calculated at po- We take m u = 5 MeV. 





Mj (MeV) 


Rj (fm) 


B 1 / 4 (MeV) 


z 


M* (MeV) 


R* (fm) 


p (CSB) 
n (CSB) 
N (SU(2)) 


937.6423 (input) 
939.6956 (input) 
939.0 (input) 


0.8 (input) 

0.8000 
0.8 (input) 


169.81 
169.81 
169.97 


3.305 
3.305 
3.295 


751.928 
753.597 
754.542 


0.7950 
0.7951 
0.7864 



CHARGE SYMMETRY BREAKING IN QMC 

Now we introduce the charge symmetry breaking in the QMC model [7,11]. The 
charge symmetry is explicitly broken at the quark level through their masses. We 
use different values for the u and d current quark masses, and the effective proton, 
M*, and neutron, M*, masses. At position f in a nucleus (the coordinate origin 
is taken at the center of the nucleus), the Dirac equations for the quarks in the 
proton or neutron bag are given by 




= 0, (7) 



where \x — r\ < R* (j specifies proton or neutron). Note that we have assumed 
that the scalar potential is common to both the u and d quarks. The nucleon 
and meson fields are calculated self-consistently by solving a set of coupled non- 
linear differential equations, derived from the effective Lagrangian density (3) with 
the proper modifications caused by the different proton and neutron (or u and d 
quark) masses in MFA. Thus, the present calculation is free from the sort of double 
counting questioned by Auerbach [12], and includes the shell effects, which were 
discussed by Cohen et al. [13]. 

Before discussing the results obtained, we again need to specify the parameters 
and inputs used in the calculation [11]. They are summarized in Table 1. The bag 
constant, B, and the z parameter are determined by the bare proton mass, after 
allowing for the electromagnetic self-energy correction +0.63 MeV, with the bag 
radius, R p = 0.8 fm, in free space. For the neutron, the procedure is the same as 
that for the proton, allowing for the electromagnetic self-energy correction, —0.13 
MeV, but using the values of B and z determined above and calculating the d 
current quark mass and the bag radius for the neutron. Thus, the u current quark 
mass (m u = 5 MeV) is the basic input parameter used to fix the model parameters 
so as to reproduce the bare proton and neutron masses in free space. We found 
171,1 — 9.2424 MeV in the present calculation. 

The coupling constants, g% and g%, are determined so as to fit the saturation 
properties of symmetric nuclear matter [11]. In Table 1, SU(2) stands for the 
parameters and inputs obtained and used for the calculation when SU(2) symmetry 
is assumed, namely m u = = 5 MeV. We then found: (g%, g£) = (5.698, 2.744) 
for CSB, and (5.685, 2.721) for SU(2). For the quark-p meson coupling constant, to 



make a realistic estimate, we here use the phenomenological value, g q p = 4.595, the 
value at zero three-momentum transfer corresponding to Hartree approximation, 
from Table 4.1 of Ref. [14]. (Note that because the QMC model does not contain 
the p-nucleon tensor coupling [8], this gives an unrealistically large value for the 
coupling constant [11].) 

Proton and Neutron Masses in Nuclear Matter 

As in Equation (5), the proton and neutron masses are again given by functions 
of a in matter, and may be expanded in terms of a at low p B 

M; = M p + (3<7«)^[2S tt/p (0) + S d/p (0)]a + 0(a 2 ), (8) 

M* = M n + (3gl)^[S u/n (0) + 2S d/n (0)]a + 0(a 2 ). (9) 

Because m u ^ rrid, the u-quark scalar charge is no longer the same as the d-quark 
scalar charge. We have therefore introduced four kinds of quark scalar charges in 
the expansion: Si/j(cr) = J v . df tp^^i/j, where % denotes u or d quark, Vj is the 
volume of j (= p or n) and ifti/j is the i quark wave function in j. Since the proton 
consists of two u quarks and one d quark, the derivative of M* with respect to 
a is given by 2S U / P + S d / P . Similarly, the derivative for the neutron is given by 

S u /n + 2S d / n . 

Taking the difference between the in-medium neutron and proton masses, we find 

K P = K - M; = A np - (3(g)[S B (0) - S p (0)]a + 0(a 2 ), (10) 

where A np = M n - M p , S n (0) = §[S u/n (0) + 2S d/n (0)\ and S p (0) = |[25 u/p (0) + 
S d / P (0)]. Here we may expect that S u /j < S d /j because the u quark is more rela- 
tivistic than the d quark in nuclear matter {m u < m d ) — note that the quark scalar 
charge is given in terms of ip q ip q in matter. Thus, we find that S n (0) > S p (0) and 
then A* p < A np in nuclear medium. 

In Figure 2 we show the neutron-proton effective mass difference calculated in 
symmetric nuclear matter, including the electromagnetic self-energy corrections. 
One notices that the mass difference becomes smaller as the density increases. 
This behavior works in the right direction to resolve the ONS anomaly. 

The ONS Anomaly in Mirror Nuclei 

Now we are in a position to show our results of the ONS anomaly in mirror 
nuclei [11]. We first present the calculated single-particle energies for 17 F and 17 
in Table 2. These mirror nuclei have a common core nucleus, 16 0, and have an 
extra valence proton for 17 F and neutron for 17 0. In order to focus on the strong 
interaction effect for the valence proton and neutron, the Dirac equations are solved 



without the Coulomb and p-meson potentials, or the electromagnetic self-energy 
corrections, and keeping only the charge symmetric a and uj mean field potentials. 
Consistently, the valence nucleon contributions are not included in the Coulomb 
and p-mean field source densities in the core nucleus. However, for the nucleons in 
the core nucleus, electromagnetic self-energy corrections and the Coulomb potential 
as well as the p mean field potential are included in addition to the a and uj mean 
field potentials in solving the Dirac equations. Results are shown for two cases 
in Table 2: calculation with charge symmetry breaking (denoted by CSB) and 
calculation performed assuming SU(2) symmetry (denoted by SU(2)). 

The SU(2) results for 17 F and 17 agree perfectly with each other as they should. 
Single-particle energies in the cores of 17 F and 17 are slightly different for CSB. 
This difference is induced by the different (effective) masses for the valence proton 
and neutron, arising from the charge and density dependence of their coupling to 
the self-consistent scalar mean field. This also causes a second order effect on the 
Coulomb and p-meson potentials through the self-consistency procedure. 

It is interesting to compare the binding energy differences between the va- 
lence proton in 17 F and neutron in 17 0. In CSB, the result gives, E(p)(ld^/2) — 
E(n)(ld^,/2) — 0.18 MeV, while the SU(2) case is zero as it should be. This amount 
already shows a magnitude similar to that of the observed binding energy differ- 




0. 70 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



FIGURE 2. Neutron-proton effective mass difference in symmetric nuclear matter with the 
electromagnetic self-energy corrections. 



TABLE 2. Calculated single-particle energies 
(in McV) for 17 F and 17 0. 





CSB 

17 F 17Q 


SU(2) 

17 F 17Q 


P lSl/2 
1P3/2 
lPl/2 


-28.800 -28.805 
-14.154 -14.158 
-12.495 -12.499 


-28.663 -28.663 
-14.032 -14.032 
-12.383 -12.383 


n lSi/2 
1P3/2 
lPl/2 


-33.367 -33.372 
-18.259 -18.263 
-16.587 -16.590 


-32.967 -32.967 
-17.918 -17.918 
-16.258 -16.258 


valence 
ld 5 /2 


p n 
-3.918 -4.099 


p n 
-3.848 -3.848 



ences. 

In Table 3, we summarize the calculated single-particle energies for the valence 
proton and neutron of several mirror nuclei (in CSB) [11]. Comparing the p- 
potential contributions for the hole states with core plus valence states, one notices 
the shell effects due to the p-potentials. These results reflect the difference in the 
shell structure, namely the hole states tend to have larger p-potential contributions 
than the core plus valence nucleon states. 

The binding energy differences obtained indicate that the prime CSB effects 
originate in the u-d current quark mass difference. The calculated binding energy 
differences give of the order of about a few hundred keV. This is precisely the order 

TABLE 3. Calculated single-particle energies of mirror nuclei. For each 
nucleus, the top row shows the single-particle energy of the valence proton 
or neutron (the orbit is also indicated). SE p stands for the contribution 
from the p-meson central and spin-orbit potentials of the core nucleus. The 
discrepancies between the experimental values and the theoretical expec- 
tations in the absence of charge symmetry violating strong interactions are 
taken from Table II of Ref. [5], by averaging over the theoretical values. 





15 0( P ) 


15 N(n) 


17 F( P ) 


17 0(n) 


lp V 2 or ld 5/2 (MeV) 


-14.397 


-14.631 


-3.918 


-4.099 


6E p (MeV) 


-0.055 


0.056 


-0.005 


0.005 


Total(MeV) 


-14.452 


-14.575 


-3.923 


-4.094 


SE = E(p) - E{n) 


5E = 


123(keV) 


5E = 


171 (keV) 




observed = 


230 (keV) 


observed = 


220(kcV) 




39 Ca(p) 


39 K(n) 


41 Sc( P ) 


41 Ca(n) 


ld 3/2 or lf 7/2 (MeV) 


-16.407 


-16.689 


-6.970 


-7.210 


SE p (MeV) 


-0.087 


0.088 


-0.006 


0.006 


Total(MeV) 


-16.494 


-16.601 


-6.976 


-7.204 


5E = E(p) - E(n) 


SE = 


108(keV) 


SE = 


228 (keV) 




observed = 


340 (keV) 


observed = 


460 (keV) 



of magnitude which is observed as the ONS anomaly [3,5]. 

SUMMARY 

Using the QMC model, we have discussed CSB in nuclear medium and calcu- 
lated the ONS anomaly in mirror nuclei, including the quark degrees of freedom 
explicitly. We stress that the present contribution to the ONS anomaly is based 
on a very simple but novel idea, namely the slight difference between the quark 
scalar densities of the u and d quarks in a bound nucleon, which stems from the 
u and d quark mass difference [7,11]. This implies that the in-medium proton-o" 
and neutron-a coupling constants differ from their values in free space and that the 
neutron-proton effective mass difference is reduced in matter. 

Our results were obtained within an explicit shell model calculation, based on 
quark degrees of freedom. They show that once CSB is set through the u and d 
current quark mass difference so as to reproduce the proton and neutron masses 
in free space, it can produces binding energy differences for the valence (excess) 
proton and neutron of mirror nuclei of a few hundred keV. The origin of this effect 
is so simple that it is natural to conclude that a sizable fraction of CSB in mirror 
nuclei arises from the density dependence of the u and d quark scalar densities in 
a bound nucleon. 

It is a fascinating challenge for the future to compare this result with the tra- 
ditional mechanism involving p — uo mixing [4] . This will involve the issue of the 
possible momentum dependence of the p — uj mixing amplitude [3,15]. In addition, 
one would need to examine whether there is any deeper connection between these 
apparently quite different sources of charge symmetry violation. 
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